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I For a large region of parameter space involving the cosmological constant and mass 

I parameters, we discuss fluctuating spacetime solutions that are effectively Minkowskian 

on large time and distance scales. Rapid, small amplitude oscillations in the scale factor 

! have a frequency determined by the size of a negative cosmological constant. A field with 

^ I modes of negative energy is required. If it is gravity that induces a coupling between the 

CN ! ghost-like and normal fields, we find that this results in stochastic rather than unstable 

, behavior. The negative energy modes may also permit the existence of Lorentz invariant 

, fluctuating solutions of finite energy density. Finally we consider higher derivative gravity 

CD ■ theories and find oscillating metric solutions in these theories without the addition of other 

^ . fields. 
-4— > , 

^ • 1 Introduction 

> : 
•1—1 , 

^ ■ One ingredient of the cosmological constant problem is the observation that the classical field 
d • equations of General Relativity do not admit a Lorentz invariant solution in the presence of a 
cosmological constant A. In this paper we will reconsider this point after relaxing the positive 
energy conditions on the energy momentum tensor. In particular we will consider the effect 
of negative energy contributions in the presence of a large A. We will suggest that it then 
becomes possible that the theory has solutions which are fluctuating on scales characterized by 
A in such a way that no preferred frame is chosen, and which appear to be Minkowskian on 
larger length and time scales. That is, the dynamics controlled by a nonzero A involves rapid 
fluctuations that are effectively averaged over by a low energy observer, and are such that no 
Lorentz symmetry breaking is induced in the low energy theory. For the moment we will simply 
comment on some consequences of such solutions, assuming they exist, keeping in mind that 
we will also be faced with issues associated with negative energy modes. 
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Although wc arc supposing that Lorentz invariance is preserved on high energy scales char- 
acterized by A, it is not explicitly preserved by nature on lower scales. It is broken by the 
energy momentum tensor T^j, of whatever system we wish to study, and in particular it is 
broken by the T^^, that defines the preferred cosmological frame. We would like to describe 
these explicit violations of global Lorentz symmetry as long wavelength perturbations to the 
high energy Lorentz invariant solution. If wc consider performing a derivative expansion of the 
original coordinate invariant theory around the Lorentz invariant symmetry point, then this low 
energy theory up to two derivatives must be General Relativity. In order for it to be consistent 
with its interpretation as an expansion about a solution that is flat on large distances, the low 
energy theory must have vanishing cosmological constant. 

We stress that given our assumptions, the vacuum energy component of T^^, does not intro- 
duce explicit Lorentz symmetry breaking, whereas the remainder of T^,^ does. The cosmological 
constant influences the dynamics of the Lorentz invariant high energy solution, whereas the 
explicit Lorentz violating effects are accounted for in the low energy theory. Nevertheless, it 
seems paradoxical that the cosmological constant, which could be considered an inflnite wave- 
length component of T^i^, has been absorbed by high energy physics. But this is a direct result 
of our assumption that A influences the high energy theory by determining the amplitude of the 
Lorentz invariant fluctuations, and thus yielding solutions that are Minkowskian in the infrared. 
This is the flip side of the statement that the low energy theory does not contain a parameter 
A which naively describes inflnite wavelengths. 

In section 2 we will illustrate the essential point with a model containing two decoupled fluids, 
with one of them having negative energy. The gravitational effects of these fluids will lead to 
a metric that is asymptotically flat in the presence of nonvanishing A. y^|A[ still introduces 
a length scale in these solutions, just as it does for dS or AdS solutions, but here it sets the 
frequency of oscillation. The scale factor in the metric oscillates, as do the curvature invariants. 
The amplitude of this oscillation can vanish for special initial conditions, but typically it does 
not. The point about oscillating solutions is that their Minkowskian appearance on larger scales 
is a universal property of this whole space of solutions. Thus these solutions have an ultraviolet 
property, their frequency, and an infrared property, their effective flatness. The linking of these 
properties is at the origin of the counter-intuitive picture we are developing. 

What is missing in these prototype solutions is Lorentz invariance, since the solutions are 
characterized by a preferred frame. To have a state that is both rapidly fluctuating and Lorentz 
invariant we need the Lorentz invariant analog of a rotationally invariant thermally fiuctuating 
system. But unlike a rotation, a boost can generate states with arbitrarily large energies, and 
thus it appears that a nontrivial finite energy Lorentz invariant state cannot exist. We will 
argue that the situation can change in the presence of the negative energy modes, since they 
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can serve to regulate the apparent infinite energies. 

Of course negative energy modes, gliosts, raise tlie more immediate issue of stability, and 
so in section 3 we study interacting ghost dynamics using massive scalar field theory where the 
masses are large compared to the A scale. A ghost instabihty can occur if the ghost couples 
to other fields, and so we explore the effect that gravity has when gravity is responsible for 
the coupling between the positive and negative energy modes. After including the effect of a 
changing scale factor on this coupling we find that solutions exhibiting unstable behavior are 
typically not present. In section 4 we consider a superposition of solutions involving different 
momentum modes and thus recover, at least formally, Lorentz invariance in the infinite mode 
limit. In section 5 we discuss a qualitative difference that occurs when a number of different 
fields are introduced, namely that the energy per mode need not be vanishing. As a classical 
model for Planck scale physics we study pure gravity theory involving higher derivatives in 
Section 6, where we find the analogs of the solutions in section 3. We conclude in section 7. 

We emphasize that in this work we are using classical analysis to model physics that is likely 
close to the Planck scale. We do this because the dynamics of interest may well manifest itself 
in classical equations. This analysis misses the quantum instability associated with a ghost PP, 
for example due to virtual graviton exchange if nothing else 0. But if the ghost mass is of 
order the Planck mass then the virtual graviton is far offshell in the deep Planckian region, and 
there is no guarantee that the perturbative quantum description holds. We assume that it does 
not. Although our discussion in section 5 will lead us quite close to making some connection to 
quantum theory, we will leave this to pursue elsewhere. 

2 A Ghost Fluid 

We start with the notion of a ghost-like perfect fluid, completely decoupled from everything 
else except gravity. Whatever the constituents of this fluid are, their coupling to gravity will 
induce interactions between the fluid and other matter. We will return to the effect of these 
couplings in the next section, but for now we simply consider the effect of a ghost fluid on the 
gravitational field. We will only keep one degree of freedom of gravity, that of the scale factor 
A(t) in the spatially flat FRW metric, ds^ = dt^ - A(t)dx^. 

It is not surprising that a ghost fluid with negative energy and pressure and a normal fluid 
with positive energy and pressure could in combination cancel a cosmological constant, for some 
fixed value of the scale factor. We will demonstrate here that this can be a stable situation 
when the cosmological constant is negative. Namely, for deviations away from this solution, 
the solutions exhibit an oscillatory behavior for A{t). The essential ingredients are a negative 
cosmological constant A < 0, a ghost fluid with equation of state parameter Wg > Q (that is 
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Pg = WgPg wlth both Pg pg nGgativc), and a normal fluid with, as we shall see, < w„ < Wg. 

Before solving the equations we describe the origin of the oscillation. If A{t) oscillates then 
it has a maximum and a minimum, and at these points the total energy vanishes since it is 
proportional to H(t)'^ {H = A/2A is the Hubble parameter). Also at these points oc —H{t), 
and thus the total pressure must be negative at the minimum and positive at the maximum. 
Starting at the minimum scale factor we have the energy + p„ + A = and the pressure 
+ P„ — A < 0; and since A < 0, /)„ dominates pg and pg dominates pn- Given that Wg > Wn, 
as the universe expands \pg\ will fall faster than p„ resulting in positive total energy. \pg\ will 
fall faster than p„ so that the total pressure becomes less negative and eventually changes sign 
as A becomes relatively more important. Both \pg\ and pn decrease relative to A, and after the 
pressure has changed sign, the total energy returns to zero, at which point the maximum scale 
factor is reached. 

The "energy" and "pressure" Einstein equations, making use of the conservation of the 
separate fluids, are 
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4 A{ty 
A{t) 1 A{tf 
A{t) 4A{t)^ 



= A-A(i)-^/'-^/'-V, + >l(i)-^/'-3/'-"r„, (1) 
= A + A(i)-3/2-3/2-«^^^^_A(i)-3/2-3/2«,„^^^^_ (2) 



Here r„ and are positive constants. We start with a trivial time independent solution where 
for convenience we set A{t) — 1. Then 

Tg = - A ( 1 + w,, ) / {Wg -Wn), 

r° = -A{l + yjg)/{Wg-Wr,). (3) 

Since the r's are positive we see again that Wg > Wn ii A. < 0. We also see that if Wg is close to 
Wn then the individual energy densities in the two fluids can be large compared to A. 
We consider a small deviation from these values of the form 

Vg _ . 3 
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1+2(1 + V 
5 = 1 + + (4) 

' n 

where £ is a small number of either sign. The following oscillating scale factor with amplitude 
l^l is then a solution, to first order in e. 

A{t) = 1 + (1 - cos{u;t))e (5) 
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The frequency is 



u; = \l~-{l + Wg){l + Wn)A, (6) 
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where A = A at this order. Here we see why A must be negative; a solution of the form (j2)) 
exists for A > and Wg < Wn, but it would not be stable under small deviations. 
At second order in e, we find that A is shifted 
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A = A 1 + — {2wn + l~Wg){2wg + l-Wn)e'] , (7) 
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1 + ^(1 + Wn)£ + ^(1 + Wn)iWn + Wg)e^, 

^ = l + \{l + wg)e + ^{l + wg){wn + wg)e\ (8) 
The oscillating A{t) picks ups a higher harmonic. 

A{t) = ! + (!- cof^{ujt))e + + Wg) [cos(u;t) - cos(2cjt)] (9) 

8 

The solution can be extended to order by adding appropriate terms to r„ and r^, so that A{t) 
picks up new terms proportional to cos(ti;t), cos(2a;t), and cos(3ci;t) respectively. The existence 
of these solutions to all orders in e can easily be checked numerically. 

For a given value of a negative A there is no fine tuning of or r„ necessary to obtain 
some solution in this class. The solutions are characterized by two properties, the averaged 
value of the scale factor and the amplitude of its oscillation, and these will adjust depending 
on the values of Vg and r„. (It was only for convenience that we expanded the solutions around 
A{t) = 1.) 

It would seem that we can simply average over these rapid fluctuations, leaving something 
that looks like Minkowski space on scales large compared to 1/uj. But the problem is that 
the effective theory below the A scale is not General Relativity. The scale factor is oscillating 
around a fixed value determined by the energy densities in the conserved fluids, and so the 
scale factor degree of freedom no longer exists in the low energy theory. In other words this 
gravitational degree of freedom has picked up a mass, given by the frequency of oscillation in 
®. 

This problem is closely related to the loss of Lorentz invariance in the low energy theory, a 
result of considering fluids with a preferred rest frame in which only spatially constant solutions 
are considered. It is also related to our assumption that the fluids are noninteracting and 
separately conserved. But in fact interactions between these fluids are induced, a topic to 
which we now turn. In section 4 we return to the issue of Lorentz invariance. 
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3 Ghost Interactions 



The previous discussion assumed that the ghost and normal fluids are completely decoupled. If 
we supppose that these fluids are composed of massive flelds then gravity will induce a coupling 
between the flelds, even at the classical level. This is because massive flelds involve rapid time 
oscillations, and these oscillations feed into the metric through the gravitational coupling, and 
this in turn induces interactions between the positive and negative energy modes. 

Before we consider this, let us flrst consider the effects of a nongravitational coupling between 
the ghost and normal flelds. This coupled system can of course suffer from a severe instability, 
where the amplitudes of both flelds can grow without bound while conserving total energy. The 
nature of the instability depends on the form of the coupling between these flelds For some 
types of couplings the instability only occurs if the amplitudes of the flelds and/or the size of 
the coupling is sufficiently large, while for other types of couplings the instability always occurs. 

The simplest example of the latter case is when a ghost scalar (pg and a normal scalar 0„ 
have a coupling 4'g(f)n, but are otherwise free flelds. We will refer to these two real flelds as the 
ghost and scalar respectively. 

C = \d,<l)ndyn - \ml(t>l - \d,(t>ad^(t>9 + ^n^g + A0„0, (10) 

Note that the mass term of the ghost is also of "wrong sign" so as to ensure oscillatory solutions. 
For solutions of the form 0„ = a{t) cos(p„ ■ x), (pg = rj{t) cos(pc, ■ x), we look at the instability 
that occurs for equal masses m„ = mg = m between momentum modes that are in resonance 
p„ = ±Pg. In terms of the resonant frequency uo = y/]?~+~m? the solution is 

a{t) = ci cos(u;+t + 6i) exp(co'_t) + C2 cos(ti;+t + ^2) exp(— cj-t) 

rjit) = Ci cos(ci;+t + 5i — K7r/2) exp(cj_t) + C2 cos(co'+t + + fi:vr/2) exp(— 

2cu! = + A2 - 

2ujI = Vlu^ + A2 + cu^ 

K = sign(A) (11) 

Thus there are exponentially increasing and decreasing solutions, and they differ in the relative 
phase of the oscillating factors in a{t) and f]{t). This feature will have implications in the next 
section. 

These solutions are very model dependent. If there are other terms in the theory that 
become important for large fleld amplitudes then the exponential growth may not survive. For 
example adding terms of the form of a kinetic energy squared, for one or both of the flelds with 
appropriate signs, can change the exponential growth into an oscillating modulation of the high 
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frequency oscillations. The amplitude of the modulation factor is inversely related to the size of 
coefficients of the higher derivative terms, and its frequency is proportional to A, at least when 
this is small. 

We now turn to a ghost+scalar system where the only coupling between the two fields is 
provided through their coupling to gravity. The two fields can only interact if A{t) is noncon- 
stant. Assuming the same form for 0„ and (pg as before, the spatial dependence drops out of 
the two scalar field equations. After performing a spatial average of the ghost and scalar contri- 
butions to Tfj_^, the two field equations along with the "pressure" Einstein equation gives three 
equations for the three quantities A{t), cr{t), ri{t). The "energy" Einstein equation remains as 
a constraint. 

Since the energy densities of the ghost and scalar have to be at least of order A, and the 
frequencies of these fields are of order m ^ rrin ~ frig, then the field amplitudes are at least of 
order ^J\A\/m. In the following we will consider this ratio to be small. Our prejudice is that A 
is not extremely small in Planck units, in which case m can typically be of order mp\. But this 
need not be the case; both ^J\A\ and m can be much smaller than mpi and our description will 
still apply. 

Our main interest is in the unequal mass case where it is clear that gravity infiuences the 
resonance condition Pg + = + fn"^ for the interacting modes. The physical 3-momenta of 
the two modes are dependent on the scale factor, so that in terms of the comoving wavevector 
of each mode we have kg/A(t) and k'^/A(t). Thus a changing scale factor, which is 

responsible for producing a coupling between ghost and scalar in the first place, can cause two 
modes {kg, kn) that were in resonance to fall out of resonance. 

Numerical analysis of the two mode case indicates that a stable, oscillating system can result 
as long as kg > kn. A{t) experiences oscillations determined by A, similar to (0). We refer to 
these as A oscillations to distinguish them from the more rapid oscillations of the scalar fields. 
The basic stability is due to fact that the resonance condition for the fields is satisfied at only 
one value of the scale factor in the A oscillation cycle. This can cause a disturbance to the 
motion, so that it is no longer strictly periodic, but the basic stability remains. The condition 
kg > kn is equivalent to the previous condition Wg > w„ in the two fiuid model. The difference 
here is that we are taking the rapid time dependence of the massive fields into account, and 
including the back reaction on the metric. The latter is evident in the pressure equation where 
the matter terms oscillate and are largely canceled by an oscillating H{t) term. The ghost 
and scalar frequencies are affected by the changing A{t) as it undergoes A oscillations, and 
this leads to a beat type pattern in the sum of pressures, which in turn feeds back into A{t). 
The A oscillations that are evident in the individual ghost and scalar energies can be quite 
nonsinusoidal. 
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Figure 1: An example of a scale factor A{t) for two pairs of positive and negative energy modes. The 
dashed lines show the values of A(t) for which a resonance occurs, = 1, rUg = .8 and A = .0001 
in Planck units. A{t) also has a small amplitude, high (order 1) frequency component, but it is not 
visible. 

4 Multi-mode solutions 

First consider a pair of modes, one positive and one negative energy, and each with a particular 
4-momentum. A Lorentz transform of such a solution will also be a solution. For an arbitrarily 
large boost the positive and negative energies can be made arbitrarily large but they will nearly 
cancel, with the remainder canceling the cosmological constant just as before. For large energies 
and momenta the masses are less significant, and the ratios of energy to pressure of the positive 
and negative modes approach a common value. This is the analog of the situation described in 
(PI) when Wn ~ Wg and where the individual energy densities were also becoming large compared 
to A. 

For these Lorentz boosted solutions the amplitudes of the individual fields remain of order 



y |A|/m. If this is small then the solutions are close to being solutions to the linearized equations. 
Thus we can expect that a superposition of the Lorentz transformed solutions would also be 
a solution, up to small corrections. We can also add pairs of modes that are in resonance at 
differing values of the scale factor. The more pairs we add, the smaller the contribution of each 
to the total that cancels A. 

If a large number of modes are excited then new pairs of modes can continually come into 
resonance as A{t) changes. But as evident in (jlip . whether a pair of modes in resonance tend to 
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increase or decrease in amplitude depends on their relative phases. If the phases of the various 
modes are effectively random, then the small changes in amplitudes caused by the short-lived 
resonances will also be random. We thus expect that the individual mode amphtudes display a 
stochastic behavior/ and by the feedback on A{t) the scale factor should also have a stochastic 
component. 

We consider the case m„ > rrig, since this ensures that the pairs of modes in resonance 
satisfy kg > kn- We have numerically studied this case with up to 4 pairs of modes with varying 
values of kg and kn- We observe the stochastic behavior along with the basic stability of the 
multi-mode system. Fig. (1) shows an example with two pairs of modes. The stability even 
survives the addition of direct interaction terms between the ghost and scalar, as long as the 
coupling is not too large. 

We can contemplate taking the limit as the number of modes goes to infinity. The phase 
space is infinite, and here we are contemplating that this whole phase space is populated. 
In other words there would be positive and negative energy seas corresponding to a smooth 
distribution of modes covering the respective mass shells. We have seen how this can be built 
up from pairs of modes with nearly cancelling energies. In the limit, the energy contributed 
from each pair goes to zero, since the total energy density is finite. That is, the energy density 
per unit phase space of these paired modes should vanish in the limit. 

The interest in the infinite mode limit lies in its possible Lorentz invariance. The point is 
to have a situation where none of the fluctuating components, scalar, ghost, or metric, displays 
a preferred frame or 4-momentum. It is this that could guarantee that the low energy theory, 
obtained after integrating out the rapid fluctuations, is described by General Relativity. 

If a collection of fluctuations is to display no preferred frame under Lorentz transformations, 
then we already know what its energy spectrum must look like. It must be proportional to the 
zero point energy spectrum of a field in a relativistic quantum field theory, where the energy per 
mode is ^huj]^. This type of classical background with all modes of the relativistic phase space 
excited, and with each mode having a random phase, has been considered in There it was 
explicitly checked that when a suitable average over the random phases is performed, Lorentz 
invariant results are obtained. In our case the amplitudes of the positive and negative energy 
modes must be such that the energy per mode is proportional to cuJJ and — cu^ respectively, 
where l<jJ, = y'k^ + mj. It is also possible for the energies to deviate from these values in some 
random way. 

The main point for us is that the sum of positive and negative energy modes can be finite. 
We are suggesting that this finite result is due to the same dynamics that we have been able 
to study numerically for a finite collection of modes. This is the dynamics that responds to a 
^The system is of course deterministic for a given set of initial conditions. 



9 



finite A to produce effectively flat space. In this way the infinite mode system is self-regulating. 
In fact it is very close in spirit to Pauli-Villars regularization, except that the negative energy 
states here are not fictitious. 

For a classical field theory without negative energy modes, the only Lorentz invariant con- 
figurations of finite energy density have constant fields. Here the negative energy modes allow 
the existence of less trivial configurations, which are the analogs of rotationally invariant states 
of fluctuating thermal systems. The obvious question then is why this unorthodox distribution 
over phase space doesn't cascade down and thermalize. The formal answer is Lorentz invari- 
ance; any thermalization implies that a preferred frame has been selected and thus a breaking 
of Lorentz symmetry. The unbroken symmetry prevents this from happening. For the same 
reason there can be no transfer of energy to or from the fluctuating background (for a related 
discussion see jSj). 



5 Many Fields 

We can analyze further how the normalization of the energy spectra is dynamically determined. 
We first consider our previous case of just one ghost field and one nonghost field with different 
masses. We introduce a momentum cutoff M and consider the approach to the M — >■ (X) limit. 
The "energy" and "pressure" components of the energy momentum tensor are 

A + y ^ {anu^l - agui) = (12) 

-A+ / - — -- [an— a„^\ = 0. (13) 

J (27r)3 \^ "3^1: '^SuiJ ^ ' 

The Oj's determine the energy per mode normalization; we have omitted k subscripts because 
we have seen that they are k independent (at least on average) due to Lorentz invariance. These 
equations constrain a„ and ag in such a way that they are vanishing as M —>■ oo. 

+1^^ + 0(M-.) (14) 



{ml - m2)M2 

= (,„g _ „.)M^ - + 0(M-% (15) 

Thus the energy per mode is vanishing in the infinite mode limit, and this result is consistent 
with our numerical findings for this case. 

But this result changes when more fields are involved. For fields we have 

A + / 7^ E«^-U = 0^ (16) 
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^lE'-.^^l = 0. (17) 




We expand the coefficients as follows 

,(2) ^(4) 



where a^-^^ and a^'^^ are linear functions of ln(M). The equations can have solutions when 

N N N 

Ear=E«SV = E«r = 0. (19) 

i=l i=l 1=1 

These three constraints are sufficient to cancel the quartic and quadratic divergences.^ Thus to 
have nonvanishing as M — > cxd there must be at least one ghost and one nonghost and > 3. 
For larger we note that if there are at least three massive fields (involving both ghosts and 
nonghosts) then the remaining massless fields can have af^^^ that are equal in absolute value. 

This shows how the energy per mode can stay finite in the M — oo limit. In fact we 
would expect that all fields of the theory should be involved, due to the interactions between 
fields. But as soon as we consider light fields, we are forced out of our classical discussion, since 
light fields are quantum mechanical. This in turn means that light ghost fields are forbidden 
since they would introduce a severe quantum mechanical instability (or there needs to be a 
very low energy Lorentz violating cutoff [21 El)- On the other hand the classical Dirac theory 
describes negative as well as positive energy modes. As far as the classical analysis of this 
paper is concerned, the ffiling of the negative energy modes of this classical Dirac theory would 
serve the same purpose as the filling of the negative energy modes of a classical ghost theory. 
The essential difference is that the Dirac theory and its negative energy sea may survive the 
transition to quantum mechanics without the introduction of instabilities. 

All of this leads to some speculations concerning the foundations of quantum field theory, 
with h, as arising from the energy per mode normalization of classical fluctuations, which in turn 
is determined by A. We are pursuing this by calculating correlation functions from random phase 
classical fluctuations on a 1 + 1 dimensional lattice, and comparing to the correlation functions 
of quantum fleld theory. We will report on these surprisingly positive results elsewhere 



^The elimination of the log divergences and the fixing of the finite pieces determines combinations of af'^ 



and a. 



(4) 
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6 Pure Gravity 



Thus far we have used fluid models and scalar field theories to model ghost dynamics. In this 
section we would like to show how similar dynamics can arise if ghosts have a purely gravitational 
origin. We again seek solutions involving classical fluctuations of the metric around Minkowski 
space, with an amplitude determined by A. Of course such solutions do not exist in pure General 
Relativity. The required new massive degrees of freedom can appear when GR is supplemented 
with higher derivative terms. Higher derivative theories are well known to have ghosts with 
Planck scale masses, which are seen when linearizing the theory j^. For a different perspective 
on these ghosts, see [Zj. 

As before, for the solutions of interest the amplitudes of the fluctuations will have to be of 
order ^|A| in Planck units, so that the energy can be of order A. We are able to flnd such 
solutions when A is small and negative. Thus the basic result for pure gravity with negative A is 
that besides the anti-de Sitter solutions, there are also solutions that are effectively Minkowskian 
at large distances. 

Our action is the general gravitational action with 4 as well as 2 derivative terms, including 
a cosmological constant. 



d^xy^ (A + kR + aR^ + bR^^R^"" 



S 



We choose to study the various physical degrees of freedom with a metric of the form 



(20) 
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~A{t) 
-B{t,z) 
-B {t,y) 





-B {t,z) 

-A{t) 

-B 





'B{t,y) 
-B {t,x) 
-A{t) 



(21) 



B{t,x) = A{t)C{t)sm{kx), etc, 

such that g^i, = rj^y + 5g^y, where dg^^ models the small amplitude fluctuations. Thus A{t) is 
of order one but it contains an oscillating component describing a massive scalar mode. C [t) is 
a graviton amplitude, containing both massless and massive degrees of freedom. The massive 
graviton is the ghost. The amplitudes of all these modes are small when A is small. The masses 
of the scalar and massive graviton modes, as determined by the terms linear in these oscillations, 
are = —1/6 and = l/(6a + 2b) respectively in Planck units. The allowed range of a and 
b is then such that the mass squares are positive, but we otherwise take the parameters to be 
order unity. The parameter k flxes the magnitude of the massive graviton 3-momentum. 
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Although the metric breaks rotational symmetry, there are remaining discrete symmetries 
under which x, y, and z are interchanged which are convenient for simplifying the field equations. 
We also simplify the equations by making use of a small amplitude expansion for C{t). A 
wave-type equation, linear in C{t), can be obtained from the lowest order terms in the off- 
diagonal space-space components of the field equations. In the other equations the leading C{t) 
dependence begins at quadratic order, and we consider the equations that result after a spatial 
averaging of these terms. This amounts to averaging the effective T^j, from the graviton modes, 
and it leaves us with the energy and pressure equations with k dependence. 

We can use the pressure equation and the graviton equation to numerically solve for A{t) and 
C{t). The energy equation should be satisfied identically during the evolution; this offers a good 
consistency check on the equations and the numerical routines. We provide initial conditions 
and then let the system evolve; if the system doesn't diverge or hit a singularity, and if this is 
true with no fine tuning of the initial conditions, then the solution can be considered classically 
stable. This is what we find for a class of oscillating solutions having a wide range of initial 
conditions — from one solution the initial conditions can be changed in quite arbitrary ways to 
obtain other solutions. A is in the range — £ < A < where e fa 0.0005. 

The rapid oscillations in A{t) and C{t) correspond to the masses of the scalar and graviton 
modes given above. The ampfitude of these rapid oscillations is proportional to ^f\A'\. There 
is also the expected longer period oscillation with a frequency of order ^JJa]. The amplitude 
of this oscillation in A{t) is strongly dependent on the initial conditions. These features and 
the resulting stabifity is similar to what we described in the scalar field model. But there the 
ghost and scalar didn't interact except through a minimal coupling to gravity, and so it is 
interesting that the interactions between the massive scalar and graviton modes that occur in 
higher derivative gravity don't disrupt the oscillating solutions. We also note that the existence 
of the massless graviton doesn't disrupt the solutions. 

The parameter k of the massive graviton appears in the equations as k'^/A{t), the square 
of the physical 3-momentum. This is constrained not to vanish. Since we have chosen the 
3-momentum of the massive scalar mode to vanish, this is equivalent to our previous condition 
Wg > Wn- We can demonstrate the role of the massive graviton by removing it altogether, 
replacing it by some perfect fiuid with equation of state p = wp, and leaving the massive scalar 
of higher derivative gravity. Now the metric is of pure FRW form. We find from numerical 
analysis that p < 0, p < 0, and w > is required, in agreement with the picture we have already 
developed. 

As another model we replace the massive graviton by a ghost-like massive scalar field with 
a nonminimal coupling ^ adjusts the coupling between the positive and negative energy 

modes. Oscillating solutions exist once again for nonvanishing ^. If we push A down to the 



13 



lower end of its range around —0.005, the amplitudes become large and the interaction between 
the modes becomes stronger. Here we find solutions having exponential expansion; these are 
the pure gravity analogs of the solutions discussed in the appendix of 

7 Conclusion 

We began by considering the effects of a nearly decoupled negative energy fluid, which in the 
presence of a negative cosmological constant and a positive energy fluid leads to prototype 
solutions in which the metric is fluctuating and is Minkowskian on large scales. We refer to 
these fluctuations as A oscillations. We also used massive scalar flelds to model the negative 
and positive energy modes; here the flelds also oscillate due to their masses. There is no flne 
tuning of the amplitude of these massive oscillations since the amplitude of the A oscillations 
will adjust appropriately. The back reaction of the mass oscillations on the metric induces 
coupling between negative and positive energy modes. But since the value of a changing scale 
factor influences which pairs of modes are coupled at any given time, in the end this does not 
induce instability for a large range of parameters and initial conditions. Instead, a stochastic 
behavior develops. 

Though the spacetime metric is Minkowskian on large scales, these prototype solutions lack 
Lorentz invariance. By considering a superposition of Lorentz transformed solutions, we can 
formally obtain a Lorentz invariant analog of a rotationally invariant thermally fluctuating 
system. The noncompact Lorentz group implies that an inflnite superposition is required, 
but the consideration of such a superposition is not new [3^. In any case we conjecture that 
the cancellation of the cosmological constant survives this superposition, as we found for the 
superposition of a few pairs of positive and negative energy modes. All these oscillating solutions 
require negative energy modes, and now the negative energy modes can be seen to be playing a 
second role in regulating the energy density of the inflnite superposition. Under the assumption 
that the inflnite superposition solution exists and is Lorentz invariant, the low energy theory 
should be General Relativity with vanishing cosmological constant. 

General Relativity is the theory of Lorentz violating perturbations to this Lorentz invariant 
state. This includes all effects of high temperatures in the early universe, phase transitions, flelds 
slowly rolling down potentials, periods of inflation, etc. But there may be several acceptable 
zero temperature ground states of the theory, corresponding for example to various local minima 
of some potential. In principle a low energy theory with an asymptotically flat spacetime could 
be set up around any of these minima, as long as the perturbations have energies less than the 
barrier energy separating the chosen vacuum from another. 

We have proceeded beyond a scalar fleld theory description of ghosts in our study of higher 
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derivative gravity theories. The higher derivative terms introduce additional Planck scale de- 
grees of freedom, one of which is a massive ghost graviton. No other fields are introduced and 
we see again the emergence of the prototype oscillating solutions. We expect that these types 
of solutions are typical of classical gravity theories with any number of higher derivatives. This 
discussion involves extremely massive negative energy modes, but in the context of light fields 
we also commented on the negative energy modes of the classical Dirac theory. 

In this paper we have tried to demonstrate some links between Lorentz invariant classical 
fluctuations, the existence of negative energy modes, and the cosmological constant problem. 
But what does any of this have to do with the otherwise successful description of nature by a 
quantum field theory? We will report elsewhere on some intriguing results that may have some 
bearing on this question. 
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